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ABSTRACT 

The  bottom  of  the  spin-wave  spectrum  is  examined  for  the  case  of  a 

thin  film  magnetized  in  its  plane.  The  bias  field  and  wave  number 

for  the  onset  of  first-order  spin-wave  instability  at  high  microwave  power  is 

calculated.  In  addition,  the  critical  thickness  d  ,  at  which  the  uniform 

c 

mode  saturation  process  changes  from  first  to  second  order,  is  obtained. 
These  results  are  in  fairly  good  agreement  with  experiment. 


Accepted  for  the  Air  Force 
Stanley  J.  Wisniewski 
Ft  Colonel,  USAF 
Chief,  Lincoln  Laboratory  Office 
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I.  INTRODUCTION 

Suhl  has  shown  that  saturation  of  the  uniform  mode  in  ferromagnetic 

resonance  will  occur,  at  a  critical  microwave  field,  via  the  nonlinear 

coupling  to  the  uniform  mode  of  pairs  of  spin  waves  of  energy  ftuW2,  and 

their  ensuing  unstable  growth.  (uu  =  Zir  x  drive  frequency.  )  However,  if 

such  spin  waves  do  not  exist,  these  "fir st-order  11  processes  cannot  occur, 

and  the  uniform  mode  does  not  saturate  until  some  higher  field  at  which 

"second-order  11  processes  appear:  growth  of  spin  waves  of  energy  h  r  .  In 

the  latter  case  at  some  value  of  the  bias  field  H,  less  than  its  resonance 

value  H  ,  spin  waves  of  energy  h U)  /2  again  become  available  and  result  in 

a  subsidiary  absorption.  Thus  the  spin-wave  energy  at  the  bottom  of  the 

spectrum  determines  the  nature  of  possible  instabilities.  As  a  measure  of 

this  bottom,  we  may  take  the  field  H  ,  defined  as  the  value  of  H  for  which 

b 

the  energy  of  the  lowest-lying  spin  wave  is  equal  to  tiw  / 2. 

kb  ° 

Calculation  of  in  a  bulk  ferromagnet  is  completely  straightforward. 
In  a  thin  film  magnetized  in  its  plane,  on  the  other  hand,  the  spin-wave 
spectrum  is  badly  distorted  at  low  wave  numbers,  essentially  because  of 
the  creation  of  non-cancelling  surface  poles.  It  is  the  purpose  of  this  note 
to  calculate  the  resultant  bottom  field  and  wave  number.  As  a  by-product, 
we  obtain  the  critical  thickness  at  which  the  transition  from  first-order  to 
second-order  processes  takes  place.  Finally,  we  compare  our  theory  with 
experimental  results. 

II.  FIELD  AND  WAVE  NUMBER  AT  THE  BOTTOM  OF  THE  SPECTRUM 


The  dispersion  relation  for  a  film  of  thickness  d  =  2E,  with  a  field  H 

2 

in  the  film  plane,  is  given  by 

=  [oj  +  0)  k^E^  -f(ju  X(kE)][m\  +  uu  k^E^  +  a1  X(kE)sin^$]  (1) 

k  h  e  m  he  m 
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where 


=  HY 

(2a) 

—  1 

X(h)  =  K  e  sinh  h 

(2d) 

2  A 

(2b) 

XU)  =  1  -  X(h) 

(2e) 

M  L 

0 

=  4ttM  Y 

(2c) 

m  o 


and  the  wave  vector  k  lies  in  the  film  plane  at  an  angle  $  to  H.  (A  -  exchange 

constant,  =  saturation  magnetization,  Y  =  (gyromagnetic  ratio)  .  )  This 

dispersion  law  is  based  on  a  "thin-film  approximation,  n  which  assumes 

excitations  that  are  constant  across  the  film  thickness.  Equation  (1)  reduces 

to  the  bulk  spectrum  as  kE  00  (x  0,  X  1 );  in  the  thin-film  limit  kE  -*  0, 

X  -*  1  and  X  kL.  The  term  to  X  results  from  surface  poles,  which  tend  to 

m~  .2 

cancel  as  kE  —  00.  The  term  uu  X  sin  $  results  from  volume  poles,  whose 

m 

density  vanishes  as  kE  0. 

The  field  H  and  wave  vector  k  at  the  bottom  of  the  spectrum  are 
b  b 

given  by  Eq.  (1  )  and 


StOr* 

k 

"aiT 


=  0 


$  =  0  , 


which  become 


to  =  [to.  +  to  h  +  to  X(k)][u\ 
be  m  b 


(!)  K  ] 

e 


(3a) 

(3b) 

(4a) 


2tu  k[u).  +  to  K  +10  X(h  )]  +  [2to  K  +  to  X7(x  )][to  + 

e  b  e  m  e  m  b 
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tO  K  ]  = 
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(4b) 
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where 


^  =  uim 
k. 

(5a) 

b 

“b = Hbv 

(5b) 

!i=kbL  • 

(5c) 

(Note  that  (d/tt  is  the  drive  frequency.  ) 

An  exact  analytical  solution  to  Eqs.  (4)  for  H  and  ul^  is  impossible  to 
obtain.  However,  good  approximations  may  be  found  in  the  thin-film  (h  «  1) 
and  thick-film  (h  »  1)  limits;  a  numerical  solution  may  be  used  for  inter¬ 
mediate  thicknesses. 


A.  Thin-Film  Limit 
With  the  definitions 


U  = 

B  = 


(!) 

m 


(1) 

e 


ID 

UD 

m 


U) 

m 


-1 

ri  =  ku 


Eqs.  (4)  become 

b  +  bX  -  B  +  [1T\  (Zb  +  X)  +  (J.  T)  =  0 
bX'  +  2r|(Zb  +  X)  +  |-iTi2(4r|  +  x')  =  0 


(6a) 

(6b) 

(6c) 

(6d) 

(7a) 

(7b) 


where  [from  the  K  «  1  expansion  of  Eq.  (2d)] 
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(8a) 


2  2  2 

X  =  1  -  UT|  +  —  q.  T|  +... 
x'  =  -l  +  j  |-iri  -  'J' r\Z  +  . . .  (8b) 


We  solve  Eqs. 
U  by  assuming 


(7)  for  b  and  r\  by  a  perturbation  expansion  in  the  parameter 


oo 


(9a) 


= )  r\\i 
Lj  n 
n=U 


n 


(9b) 


Then,  with  the  help  of  Eqs.  (6d)  and  (8),  we  equate  to  zero  coefficients  of 
each  power  of  |a  in  Eqs.  (7)  to  obtain 

b  =  7  (  V 1  +  482  -  \  )  T)  =  T  (  1  -  /—  ■  =r) 

°  2  ^  J  °  4  ^  Vl  +4BZj 

b  =  T)2  r\  =  ~  r}2  (1  -  18  r|  +  24  T]2)  (10) 

1  o  1  3  o  o  o 

b  =  T)  T)  etc.  , 

Z  O  1 

to  any  desired  order. 

In  general,  the  condition  for  validity  of  this  solution  is 


|U  «  1 

l  •  e  •  , 

L  «  r  —  -  =  X  =  50  k* 

Z  o 

V  ZttM 

o 


For  Permalloy  (A  =  l0_Cl  erg/cm;  4ttMq  =  10  oe). 


(1  la) 


(lib) 
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However,  if 


4£Z  «  1 

(JU 

—  «  =-  3  0  kMcps 

T T  ZTT 


(12a) 


(12b) 


then  from  Eqs.  (6)  and  (10)  we  obtain 


a’b  = 


UJ 


U) 


[1  -  B2  +  2B4  + 


+  — a(l-6B  +...)  +  —  o  +  .  .  .  ]  (13a) 
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H 


3BZ  + 


+  6a  + 


where 


M.6 


Z 

U) 


ai  a' 
e  m 


The  conditions  for  validity  are  now  (IZb)  and 

a  «  Z 


i.  e.  , 

L  «  —  4JVA  ^  500  A  + 

(« 


(13b) 


(14) 


(15b) 


(15b) 


The  line  at  the  left  in  Fig.  1  shows  H^,  as  given  by  Eq.  (13a),  as  a  function 
of  thickness  for  Permalloy  with  uu/’Y  =  1530  oe.  The  corresponding  values  of 
k^  =  K  /  L  are  plotted  in  Fig.  Z . 

B.  Thick-Film  Limit 


For  this  case,  we  let 


6 

C 


For  Permalloy  at  w/n  =  8.  5  kMcps. 


(16a) 

(16b) 


5 


(»o)  qH  woiioa  wnyiOBdS  acj  cmid 
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Fig.  1.  Field  for  bottom  of  the  spin-wave  spectrum  versus  film 
thickness.  Theoretical  curves  are  based  on  the  nominal  values 
A  =  1  0“6  erg/cm,  4ttMg  =  1  0^  oe ,  uu/V  =  1  53 0  oe .  The  experiments 
were  performed  on  83%  Ni-17%  Fe  films  at  about  8.  56  kMcps. 
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thickness  at  H  =  H  Film  parameters  are  the  same  as  in 
Fig.  1.  b 


so  that  Eqs.  (4)  become 

t,  (b  +  bX  -  8  )  +  6£  (2b  +  X)  +  6  =  0 
bx'c3  +  &X'G  +  26 2C2  (2b  +  X)  +  463  =  0 


where  ,from  the  r  »  1  expansion  of  Eq.  (2d)] 


x  “  r6c 


X 


1  2  2 

-Lf>  Q 
2  ^ 


Proceeding  as  before,  we  assume  solutions  of  the  form 


and  find 


or 


■i 


b(lV 


n=  0 


G  =£  G (n)6n 

n=0 


b(0)  =  B 

Jl) 


(2) 


3_ 

4 

1 


88 

(3)  1 


486 


c(0,  =  ; 

36 


G 


(2) 


1 


66 


e  tc . 


,,3  1  Z  13 

“'b  =  U,(l  '  4  C  +  ?C  +  4?P  + 


1  1  ^  1  2 
H=  —  (1--TP+  —  P  +  ...) 


Z  6 


1Z 


where 


6  1  ,  Z 

P  =  =  —  (a1  a'  ) 

p  a,’  e  m 


(17a) 

(17b) 

(18a) 

(18b) 

(19a) 

(19b) 


(ZO) 


(2  1a) 
(21b) 

(2  Z ) 
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The  conditions  for  validity  of  this  solution  are 


i-  e.  , 


and 


p  «  1 


L  »  (  4ttVZM  A  -  800  i\ 

Va'/  o 


6«r 


i.  e.  , 


L  »  zji  \  =“  1  50  A 


(23a) 

(23b) 

(24a) 

(24b) 


Note  in  the  bulk  limit  L  —  00  that  p  0  so  that  0).  uu  and  k  0.  H  ,  as 

b  b  b 

given  by  three  terms  of  Eq.  (21a),  is  shown  by  the  line  on  the  right  in  Fig.  1  , 
and  k^  as  determined  from  two  terms  of  Eq.  (21b)  is  shown  in  Fig.  2. 

C .  Intermediate  Case 

To  find  numerical  solutions,  it  is  convenient  to  solve  Eq.  (4b)  for 
2 

UU  +  uu  K  and  substitute  this  into  Eq.  (4a),  to  obtain 
b  e 

(3 ^  (piX*  +  4kT  +  2  h (pX'  +  2  h  )  =  0  (2  5a) 


U) 


b 


UU  H 

e 


(  2uX 
V  piXr  +  4  K 


(25b) 


Eq.  (25a)  may  then  be  solved  graphically  for  k.  H^,  as  determined  from 
Eq.  (25b)  is  shown  by  the  dashed  line  in  Fig.  1.  (Note  that  the  interpolation 
between  the  two  approximate  solutions  is  smooth.  )  In  Fig.  2  the  exact 
numerical  result  for  the  wave  number  k^  is  shown  by  a  dashed  line. 
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III.  CRITICAL  THICKNESS 


The  transition  from  second  order  to  first  order  instability  occurs  at 

a  critical  thickness  d  =  ZL  such  that  U)t*  at  the  bottom  of  the  spin-wave 

c  c  k 

spectrum  is  equal  to  half  the  uniform  precession  frequency,  or 


a,2  =  (  2C\Z  H  (H  +  4tt  M  )  . 

V  2J  o  o  o 

(2  6a) 

H  =  H  (d  )  . 
o  be 

(26b) 

The  critical  thickness  is  conveniently  found  from  a  plot  of  H  vs.  d  (see 

b 

Fig.  1).  However,  an  approximate  analytic  expression  for  d  may  be  ob¬ 
tained  from  Eqs.  (26)  and  (2  1a).  Using  an  iteration  method,  we  find 


d  =  26  \ 

c  c  o 


(27a) 


where 


6  =  i  \  \  +  - 

c  3  (1  +  2p)  l  9 


(l-3g)(2  +  133)  f 

(1  +  23)Z 


j  +  I  ( 1  - 3 P ) (2  +  133)  2 

(1  +  26)2  J 


+  J_  d-3e)(i-7e) 

81  (1  +  2g)3 


•} 


(27b) 


provided  3  <  1/3.  If  0  >  1/3,  first-order  processes  cannot  occur  for  any 

3/2 

thickness,  including  bulk  samples.  In  Fig.  3  we  plot  6  ,  as  given  by 

c 

three  terms  of  Eq.  (26b),  vs.  0 ,  with  scales  at  top  and  right  for  Permalloy. 
An  interesting  feature  is  the  existence  of  an  absolute  minimum  critical 
thickness  given  by 

(d  )  .  =  61.  6  \  =  3080  k'".  (28) 

c  mm  o 
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FREQUENCY  (kMcps) 


0  2  4  6  8  10  12  14  16  18  20 


Fig.  3.  Normalized  inverse  critical  thickness  versus  normalized 
drive  frequency.  The  scales  at  right  and  top  are  computed  for 
A  =  10"6  erg/cm,  4ttM0  =  10^  oe. 


CRITICAL  THICKNESS  dr  (cm  x  to 


IV.  COMPARISON  WITH  EXPERIMENT 


Comly  and  Penney  have  measured  at  8.  56  kMcps  in  Permalloy 

e  o 

films  varying  in  thickness  from  150  A  to  7400  A.  Their  results  are  shown 

in  Fig.  1,  and  are  in  fairly  good  accord  with  theory.  The  discrepancy 

probably  results  from  the  approximations  used  to  derive  Eq.  (1):  the  true 

normal  modes  will  vary  somewhat  across  the  film  thickness,  lowering  the 

spectrum  and  therefore  raising  H  ,  particularly  for  the  thicker  samples. 

b 

e 

Comly  and  Penney  also  found  2050  <  d  <  3250  A;  the  theoretical  value  of 

c 

3120  A  is  consistent  with  these  limits.  A  more  complete  account  of  these 

4 

experimental  results  will  be  published  in  the  near  future. 
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